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The traditional magnetic storage mechanisms (both ana- 
log and digital) apply an external field signal H(t) to a hys- 
teretic magnetic material, and read the remanent magneti- 
zation M(t), which is (roughly) proportional to H(t). We 
propose a new analog method of recovering the signal from 
the magnetic material, making use of the shape of the hys- 
teresis loop M(H). The field H, "stored" in a region with 
N domains or particles, can be recovered with fluctuations 
of order 1/N using the new method — much superior to the 
l/V~N fluctuations in traditional analog storage. 

PACS numbers: 85.70.Li, 75.60.Ej, 85.70. Ay, 75.60.-d 



A. Introduction 

How can one best store a song on a magnetic tape? 
The traditional analog method converts the sound signal 
into a magnetic field H{t), and then uses it to magne- 
tize the tape, being pulled at a velocity v. The remanent 
magnetization M (the magnetization left on the tape af- 
ter the field has dropped to zero) is roughly linear in H: 

M(x) = M(vt) = C x H{t) + C 2 H 2 (t) + n(x). (1) 

Here Ci represents the nonlinearity of the remanent mag- 
netization at high fields (when recording, one turns down 
the gain until the needle during the loudest sections stops 
moving into the red), and rj represents noise (one turns 
up the gain as far as possible so quiet portions don't hiss). 
(Actually, there are nonlinearities between the remanent 
magnetization and the signal H(i) for low fields as well 
(not shown in equation (H)). When real magnetic tapes 
are recorded, the signal H(t) is convolved with a high- 
frequency, large amplitude signal [Q-Q in order to remove 
these distortions (ac-biasing). This however does not af- 
fect the new method for analog storage that we propose.) 

Two other excellent methods have been developed to 
cope with the noise and nonlinearity in the remanent 
magnetization. (There are other sources of noise in a 
magnetic recording and reading process, e.g., interfer- 
ence, electronic noise, and head noise j||]^J. In this pa- 
per, we will address only the noise relevant to the mag- 
netic material: the magnetic noise.) Dolby noise reduc- 
tion does a nonlinear transformation of H(t) to boost 
quiet sections and dampen loud sections: the inverse 
transformation is applied at playback. Digital record- 
ings are even more effective. The signal can be (linearly 



and accurately) encoded as a stream of bits, and these 
bits can be recorded and reproduced without noise or 
distortion. 

How does the analog method compare to digital record- 
ing, in terms of the amount of information that one can 
store on a given piece of magnetic tape? One important 
|l],|]J(| source of noise is the lumpiness of the irreversable 
magnetization changes in the material. Magnetic tapes 
are often made up of single-domain particles; if the field 
is strong enough, some particles will rotate their magne- 
tizations to the crystallographic axis closest to the direc- 
tion of the field p] . Other materials with large domains 
will magnetize through the depinning of sections of their 
domain walls, which (roughly) jump from one pinning 
center to another. For our purposes, these details aren't 
crucial — we mostly care about the number N of these 
pinning centers or particles. We will refer to these lumps, 
imprecisely, as domains. 

Having more, smaller domains leads to a higher in- 
formation density. Averaging over N domains will re- 
duce the fluctuation in the average magnetization by a 
factor of 1/ yN (presuming the interactions between do- 
mains are not important). Thus, the number of different 
values of H that can be distinguished in the remanent 
magnetization scales like y/~N. If we subdivide the slice 
into Q portions, and magnetize each portion separately, 

we can store a/ N/Q^ different signals. Optimizing, we 
find Q = N/e, and we store e N / 2e distinguishable sig- 
nals. This is precisely what makes the binary "digital" 
recording so effective: 2® strings of Q bits are stored, and 
our formula suggests that four domains can store one bit 
(Q = N/4). Of course, substantial error correction would 
be needed in order to keep the accuracy high at this scale! 

There are times when one is stuck with an analog sig- 
nal. Important recordings have been made with these 
outdated methods (Beatles' masters) and potentially one 
might want to reconstruct signals imposed by natural 
processes (reconstructing the stress history of a plasti- 
cally deformed material). We show here that one can 
do substantially better than the traditional analog re- 
trieval, by using the portion of the magnetization curve 
M(H) near the applied field H S i gna i (e.g., by the tape 
head during recording), rather than just the remanent 
magnetization at zero field M(0). We discuss the advan- 
tages within the context of two models: the Preisach |^J|] 
model of non-interacting hysteretic domains (which de- 
spite its simplistic assumptions is a standard tool fllCfl in 
the engineering community), and the zero-temperature 
random-field Ising model (RFIM) IOJl2|] (a more real- 



1 



istic model incorporating nearest-neighbor couplings be- 
tween domains with different threshold fields). In the 
end, we improve our resolution by y/~N (from the l/>/~N 
resolution of the remanent magnetization to 1/N), which 
for a typical slice of magnetic tape with N — 10 6 do- 
mains per wavelength (see conclusion) produces a large 
improvement in fidelity. Our method also suffers much 
less from nonlinearity: although the fidelity decreases at 
large magnetizations, the signature tracks the applied 
field directly. There is a drawback to our new method 
though: by measuring the response curve, the original 
signal is necessarily erased. 

We should mention another method for dealing with 
the random noise in magnetic films that has been devel- 
oped recently jl3| by Des Mapps of Plymouth Univer- 
sity in England. Instead of the usual two heads used 
in recording (one for demagnetizing the magnetic ma- 
terial and the other to record the signal), a third head 
is added, which reads the signal right after it has been 
recorded, and sends it to a computer for analysis. Since 
the computer "knows" what the initial signal was, it can 
adjust for the inherent noise in the magnetic material, 
record the now "modulated" signal, and leave a signa- 
ture of what it has done. This provides the information 
to the "reading" head of how to compensate, during the 
reading process. Similar techniques have been indepen- 
dently developed by Indeck and Muller from Washington 
University in St. Louis Q. 



(or minor hysteresis loops) (see fig. H) 
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FIG. 1. Magnetization M of a ferromagnet as a function 
of the external field H . The hysteresis curve appears due to a 
lag between the magnetization and the field. Ms denotes the 
saturation magnetization, Mr the remanent magnetization, 
and He the coercivity (see text). The subloop a — > b — » e — > a 
shows the return-point memory. The system comes back to 
the same state a, as the external field is switched off and then 
back on. 



B. Hysteresis Loops, Subloops, Kinks, and 
Return— Point Memory 

We review briefly the various kinds of M(H) hystere- 
sis curves relevant to our discussion. (There are many 
different curves, of course, since the response depends 
on the magnetic history of the material.) A ferromag- 
netic material has the property that its magnetization M 
lags behind the external magnetic field H, as the field is 
changed (fig. |l|) . This is called hysteresis (which means to 
lag or fall behind). The largest magnetization the mate- 
rial can have (by aligning all the magnetic domains in the 
direction of the external field H) is called the saturation 
magnetization Ms- When the field is switched off, the 
remaining magnetization is the remanent magnetization 
Mr, while the field necessary to bring the magnetization 
to zero is called the coercivity He- Variations in the 
values of these properties in different ferromagnets make 
magnetic materials useful for different applications. For 
example, magnetic recording materials have high rema- 
nence and coercivity to prevent unwanted demagnetiza- 
tion 0. Therefore, magnetic materials used in recording 
will in general have "square" hysteresis loops. By sweep- 
ing from very small to very large magnetic fields, one 
explores this saturated, "outer" hysteresis loop M^ uter 
(also called the major hysteresis loop): any other field 
history will typically be discussed in terms of subloops 



Before recording, the magnetic tape is demagnetized. 
This involves imposing a slowly decaying, oscillatory field 
H (t) which leaves the material in a well-defined, repro- 
ducible state with no remanent magnetization (up to the 
noise). Analog recording takes us from this particular 
demagnetized state to a magnetized state under an exter- 
nal field H signal- in real recording, this is done by adding 
an oscillating field to the signal (see introduction), but 
initially we will consider a monotonically increasing field, 
leading to an increasing magnetization Mi mpr i nt (H) . Re- 
leasing the external field, the magnetization M re i ax again 
drops, but to a non-zero remanent M R l9na (figure |^). 
Measuring ]\/[^9 nal gives information about H S i gna i . 

We will also be interested in M measure (H) curve, 
formed by starting from the magnetized state M R lgnal 
and raising the field again. As one sees from figure |^, 
as the external field is decreased from the original sig- 
nal H s ig na i and then increased again, the magnetization 
forms a subloop. As the external field passes H S i gna i, 
there will be very generally a kink in the M measU re(H) 
curve (figure ||). This follows in a direct way, for exam- 
ple, for models exhibiting return point memory |l4|| (also 
known as wiping out B). 
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FIG. 2. The external field is increased to H signa i and then 
decreased to zero (dashed curve). In the traditional analog 
storage, the information is stored in the remanent magnetiza- 
tion M% anal . In the new method, which uses the return-point 
memory property, the information is stored in the field H signa i 
itself, read by increasing the external field from zero and find- 
ing a "kink" in the M-H curve (solid curve). 

For these systems, the subloop closes exactly: the sys- 
tem returns at H signa i to precisely the same state it 
was in at the peak of its recording field (wiping out all 
information about the excursion to lower fields). The 
curve M measure (H) above H signal thus necessarily ex- 
tends smoothly the original curve Mi mpr int{H), while 
below H S ig na i it disagrees with Mi mpr i nt : hence it must 
have a non-analyticity at H S i gn ai- Magnetic recording 
materials often wipe out rather well: decreasing the field 
from H S ig na i (and hence repeating the loop) will repro- 
duce the same subloop (including even the noise) to a 
good extent [ p^[ . On the other hand, magnetic materi- 
als — especially those like spin-glasses with important 
antiferromagnetic couplings — can exhibit "reptation" , 
where repeated cycles lead to a slow drift in magnetiza- 
tion [jl7) . Many other systems (e.g . martensites prepared 
with parallel twin boundaries helium capillary con- 
densing in a porous material H|, and superconductors 
in external magnetic fields [gO ) can exhibit the return- 
point memory to various extents; see the work by Amen- 
gual et al. |L8| for an experiment reproducing incredi- 
ble fine structure within repeated loops. Reference Jll]] 
discusses three conditions Jl4[ (partial ordering, no pass- 
ing, and adiabaticity) which suffice to produce a perfect 
return-point memory. The models studied here possess 
these properties: antiferromagnetic couplings violate "no 
passing" . 

For the purposes of this paper, the kink in 
M me asured(H) is easily explained. As the field is raised 
a second time, the domains which flipped upon the first 
rise during imprinting and which did not flip back upon 
relaxing are not active. When H crosses H S i gna i, new, 



"virgin" domains are explored: more domains will flip 
per unit field, leading to a discontinuity in the slope. It 
is precisely this slope discontinuity in the magnetization 
curve which we propose to use in information storage. 
Since the slope discontinuity is at H S i g nai, we are saved 
from the nonlinearity and distortion which occurs with 
only measuring the remanent magnetization. We will see 
in the next section that we also suppress the noise greatly. 

C. The Preisach and Random Field Ising Models 

The comparison between the "traditional" magnetic 
analog storage method, in which the magnetization car- 
ries the information, with the new method where the 
whole hysteresis loop is retrieved, and the information is 
stored directly in the magnetic field, is done using the 
Preisach and random field Ising models. The Preisach 
model J8|,^| consists of a system of independent domains, 
each with an upper field H u and a lower field Hd, at 
which the domain changes sign (or direction). Each do- 
main therefore has its own square hysteresis loop (see 
figure ||), and the superposition of all these loops gives a 
hysteresis loop as in figure [j]. In general, the upper and 
lower fields are not equal and opposite, but the magneti- 
zation per domain mo is assumed constant. The Preisach 
model has the necessary properties for its hysteresis loop 
to exhibit return-point memory and is quite useful be- 
cause of its simplicity, but since there are no interactions 
between the domains, the modeling of real systems is 
more limited. The random field Ising model |rj]]l^1 , on 
the other hand, includes nearest neighbor interactions. 
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FIG. 3. Hysteresis loop for one domain in the Preisach 
model. The field at which the domain flips up is H u and 
the field at which it flips down is Hd- A superposition of 
many such domain hysteresis gives a hysteresis curve as in 
figure [j]. 



We first review the elements of the Preisach model 
that we will need. To calculate the magnetization in the 
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Preisach model, a weight function p(Hd, H u ) (H u > Hd) 
needs to be denned, p has the units of magnetization per 
unit field squared. For systems with time-reversal sym- 
metry, p is symmetric around the line H u — Hd in the 
(Hd, H u ) plane [|]|)| (see figure ||). If all the domains are 
pointing down, and we increase the external field from a 
large negative value to H S i gna i in figure the domains 
whose field H u fall in the shaded area of figure | will flip 
up. The magnetization of the system is then given by: 

M = J 'J p(H d , H u ) dH u dH d - J 'J p(H d , H u ) dH u dH d 

(2) 

where the first integration is over the area where the do- 
mains are pointing up (U), and the second is over the area 
where the domains are pointing down (D). If the system 
follows the path depicted in figure |l| (from a large neg- 
ative field to H signal to Hb to H e ), we obtain a "step" 
of flipped spins in the (Hd,H u ) plane (figures ||(a,b)). 
(Several subloops would give a "staircase" .) The magne- 
tization is again obtained using equation 
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FIG. 4. The triangular region is the region in (Hd,H u ) 
space where the weighting function p is defined (outside of it, 
p — 0). The shape of this region is in practice much more 
general, and depends on the properties of the magnetic ma- 
terial tell, but the region is usually symmetric around the 
H u = —Hd line. As the field is raised from a large negative 
value of the field to H s i gna i (see fig. the domains in the 
shaded region flip up. The magnetization is obtained from 
equation (b|). 
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FIG. 5. (a) The (Hd,H u ) plane after the field has been 
increased to H signa i (from a large negative value), and then 
decreased to a value Hb (see fig. ^). The shaded area corre- 
sponds to spins (domains) that have flipped (and not flipped 
back), (b) The field has now been increased to a value H e . 
More spins have flipped, and we find a "step" in the Preisach 
plane. Several subloops would give a "staircase" (see for ex- 
ample figure |b). 
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FIG. 6. (a) The ac magnetic field used to demagnetize the 
tape before storing the analog signal, (b) The demagnetiza- 
tion process in the (Hd, H u ) plane. The "staircase" occurs 
after taking the system through smaller and smaller loops. 



The zero-temperature random field Ising model has 
the following Hamiltonian {n||^] : 

H = - JijSiSj - + hi)si (3) 

<i,j> i 

where is the nearest neighbor interaction between 
spins (domains) Si and Sj (we set all = J = 1), H is 
the uniform external magnetic field, and hi is a random 
field at site Si given by a Gaussian probability distribu- 
tion. The dynamics is such that a spin Si will flip when 
its "effective" field = J £^ ■ Sj+H + hi changes sign. 
A spin that flips can trigger other spins to flip due to 
the interaction between nearest neighbors: avalanches of 
spins are possible. This model gives rise to a hysteresis 
loop and has the return-point memory property . 



We use the Preisach model with a constant magnetiza- 
tion m-o per domain, to calculate the relative fluctuations 
in the signal for the "traditional" analog magnetic stor- 
age and the new method. The weight function p can 
be written as N mo p with p being the probability distri- 
bution for the domains in the (Hd, H u ) plane. We also 
simulate magnetic systems of different sizes using the ran- 
dom field Ising model. These are then compared to the 
analytical results. 

D. Traditional Analog Magnetic Storage 

In the traditional analog recording, the tape (system) 
is first demagnetized [^lj by applying a strong ac field 
which is gradually reduced to zero (figure |^a) . In figure 
||b, the demagnetization process is show in the (Hd, H u ) 
plane for the Preisach model. In the limit of a very fine 
"staircase" (the ac field needs to drop off very slowly), 
the magnetization of the system become zero. 
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FIG. 7. Recording the signal as the remanent magnetiza- 
tion Mft 9na , starting from a demagnetized system (in the 
Preisach Plane). The shaded area represents the system at 
M = Mft 3nal and H = 0. The average value of the remanent 
magnetization j\/j^9 nal can b e calculated using equation (0). 



The tape is now ready to be recorded. The exter- 
nal field is raised from zero until a value H S i gna i (below 
the saturation value of the system), and then switched 
off. At zero field, the tape stays magnetized with the 
remanent magnetization } l /[ s ^9 nal (fig ur e |^). The mag- 
netization j\/f^9 nal can De calculated from figure using 
equation (||). However, we can calculate the magneti- 
zation differently if we notice that the probability for a 
domain to be pointing up is given by: 
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p = JJ^p(H d ,H u )dH u dH d (4) 

where the integration is over the shaded region in figure pi. 
For N independent domains, the probability P(n;N,p) 
for observing n up domains out of a total of N is given 
by the binomial distribution p2| : 

P(n;N,p) = *p-— p «(X- p )»-» (5) 
n! (iv — n)\ 

Then, the average number of domains that are up is Np, 
with a root mean square (rms) of y i\Tp (1 — p). 

Since we are interested in the measurement of the mag- 
netization, the average magnetization will be (2Np — 
N) m,Q. If we take the rms to be our measure of the fluc- 
tuation around the average, then the fluctuation relative 
to the saturation magnetization is: 

\AM% gnal \ = mo2y/Np(l- P ) = 2 ^pJT^pj 
\M S \ " Nm y/N U 

where p has a value between zero and one. The size 
1/vJV of this relative fluctuation limits the amount of 
information that can be stored using the traditional ana- 
log magnetic storage, in a system with TV domains. 

E. New Method for Analog Magnetic Storage 




FIG. 8. Storage of information in the field H a i gna i using 
the return-point memory effect, as seen in the Preisach plane. 
The shaded area represents the system at M = M^ 3nal and 
H = 0. As the field is increased, spins in the cross-hatched 
region flip up. At H a i gna i there is a discontinuity in the num- 
ber of spins flipping per unit field. This appears as a kink in 
the M-H curve. The integral of the weight function p over 
the black horizontal strip is (/' 2 ' — f )> an d over the white 
strip is f^ 1 '. 



With the new method for analog storage, the informa- 
tion is stored in the value of the external field H S i gna i at 
which the kink in the hysteresis occurs. Similarly to the 
traditional analog storage, we have a limitation on how 
well the information can be retrieved given by a value 
proportional to AH stgna i/H c , where AH signa i gives the 
fluctuation around the value H S i gna i, and He is the co- 
ercivity. 

As before, we start with the tape (system) demag- 
netized, and increase the field up to a value H B i gna i, 
smaller than the saturation field. The external field is 
then switched off (fig. ^). The system will be magnetized 
with a remanent magnetization M^ . When the in- 
formation needs to be recovered, instead of "reading off" 
the value of the remanent magnetization ]\{^9 nal ) ^ ne ex _ 
ternal field is increased until a kink in the M-H curve is 
found. If the system exhibits the return-point memory 
property, the field at which the kink occurs is H S i gna i . In 
the (Hd, H u ) plane (figure ||), the kink is seen as a dis- 
continuous increase in the number of domains flipped as 
the field is increased past H S i gna i . (If the original signal 
is ac-biased, as in audio signals, there is a discontinuity 
in the Preisach plane Q at a field shifted up from H S i gna i 
by about the amplitude of the ac-bias (see figures |](a,b). 
The retrieval of this value is otherwise analogous to the 
presentation that follows.) 



To observe the kink, we can take the derivative of the 
hysteresis curve with respect to the field H and observe 
the discontinuity. In general the data will need to be 
smoothed over some range AH , which will help in finding 
the discontinuity in the slope, but will introduce fluctua- 
tions around H S i gna i of the order of AH. The discontinu- 
ity in the slope can be observed if the difference between 
the number of domains that flip in a range AH, below 
and above H S i gna i , is larger than the standard deviation 
in the number of domains flipped in AH above H S i gna i : 

An\ 2) - An\ 1] > ^Np( 2 )(l-p( 2 )) (7) 

where the superscripts (1) and (2) indicate measurements 
below and above H S i gna i respectively, and p^ is the 
probability that a domain flips from down to up in a 
range AH, just above H S i gna i. (In general we should re- 
quire that the difference be larger than the fluctuations 
just above and below H S i gna i, but since the fluctuations 
below H s i gna i are smaller, we can use equation (^).) Note 
that Np^ = A A^^ 2 -* . If the interval AH is small enough, 
that the slope / = dM/dH of the M-H curve can be 
considered close to constant, we have: 

(f(2) _/■(!)) fm ~ 

— J - — >-AH>k J - — AH(l-p( 2 )) (8) 

m y m 

Thus the uncertainty in the measured H S i gna i is: 
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FIG. 9. (a) Magnetic field as "seen" by a magnetic tape 
moving past a recording head with an ac-biased field super- 
imposed on a constant field H signal- (b) The Preisach plane 
after the ac-biased signal in (a) has been "stored" . If the field 
is now increased (from Hu — 0) , there is a discontinuity in the 
number of domain flips per unit field as we pass the "large" 
step. That value of the field corresponds to approximately 
H S ig n ai shifted by the amplitude of the ac-bias. 

From figure ||, in the {H d , H u ) plane, is 
Nm Q f_H^ at P{ H d, H S i gna i) dH d and the difference 
(/ (2)_ / (i )) isNmo /° ^ p(H d ,H signal ) dH d . Then, 



the fluctuation in the field H S i gna i relative to the coer- 
civity H c is: 

\AH stgnal \ 1 (1 -P {2) ) I H hZL P(H d ,H slgnal ) dH d 
lHcl N \H c \(£ H3ignal p(H d ,H signal )dH d y 

(10) 

The ratio multiplying 1/N in equation ( |lC)| ) is of order 
one as long as the signal is not too small. For small 
H signal, this ratio diverges since in the Preisach model 
near (M = 0, H = 0), the M-H curve is quadratic, and 
the difference (/ ( - 2 - ) — f^) between the two slopes is neg- 
ligible. This divergence can be avoided if the signal is 
stored after the system has been saturated (instead of 
starting with a demagnetized system). Therefore, away 
from (M = 0, H — 0), the "number" of fields that can be 
used to store information scales as N. 



F. Random field Ising Model Simulation Results 

In the previous two sections, we have obtained the scal- 
ing with the system size N of the relative fluctuations in 
the magnetization (for the traditional storage method) 
and the field (for the new storage method). The analy- 
sis was done for independent domains (spins). We now 
simulate the storing and reading process for both meth- 
ods, using the random field Ising model, which includes 
nearest neighbor interactions. 

For the traditional storage method, we increase the 
field up to a value H S i gna i (from a large negative value) 
and then turn it off. We then measure the average mag- 
netization AI(H — 0) for up to one hundred initial ran- 
dom field configurations, and measure the standard de- 
viation AMg 9 "" 1 . We define the relative fluctuation as 

AM R i9na /Mri where Mr is the remanent magnetiza- 
tion, and is equal to about 0.92 (within 4%) for a disor- 
der of R — 3 (recall that R is the standard deviation of 
the Gaussian distribution of random fields). 

With the new method, we store the information the 
same way, but instead of reading off M(H = 0), we in- 
crease the external field until a kink in the magnetization 
curve is found. The field at which the kink occurs should 
be the field H S i gn ai- Figure |lO| a shows the reading pro- 
cess. We define the relative fluctuation as the difference 
between the field read off, H read , and the "real" field, 
H signal, divided by the coercivity He (which is about 
1.21 for R = 3). 



7 



0.0 




-0.8 1 ■ 1 ■ 1 ■ 1 

0.0 0.5 1.0 1.5 

Magnetic Field (H) (a) 




0.98 1.00 1.02 

Magnetic Field (H) (b) 

FIG. 10. (a) Simulation of the reading process for the new 
method of analog storage, for a 100 3 system size and disorder 
R = 3. The kink in the M-H curve corresponds to the field 
H signal ■ (b) Number of spins nipping at the field H (avalanche 
size) for the data in (a) near H signa i (which is one). The field 
H a i gn ai is found when a threshold avalanche size (here 13) is 
reached. 

To find H reac i, we note that the smaller slope of 
dM/dH inside the subloop reflects the smaller size of the 
jumps in the magnetization, or avalanches. Therefore the 
field H rea d is the field at which some "threshold" mag- 
netization jump (or avalanche size) is reached (fig. [To[b). 

Figure |ll| shows the results of our simulation. The 
diamonds correspond to the relative fluctuations in the 
field as defined above, with a threshold of 13 spins in 
an avalanche. Note that the behavior follows the 1/N 
scaling (solid line), while the relative fluctuations in the 
magnetization (squares) follow the 1 / y/~N scaling (dashed 
line). The simulation was done for 20 3 , 30 3 , 50 3 , 80 3 , and 
100 3 spins. The figure suggests a crossover at a system 



size of 100 spins, below which the relative fluctuations 
for the magnetization will become smaller than for the 
field. 



























































O AH.,„,,/H C 








□ AM R S """"/M R 




i 




o 



System size 

FIG. 11. Relative fluctuations for the field (diamonds) 
and the magnetization (square), at several system sizes. The 
solid line shows a 1/N behavior, and the dashed line a TV -1 / 2 
behavior, where N is the system size. 

G. Summary and Conclusion 

We have shown, using the Preisach model of non- 
interacting domains, that the new method of analog stor- 
age which uses the property of return-point memory, 
gives fluctuations in the signal that are smaller than the 
ones found for the traditional method. The difference is 
approximately given by a factor of vN, where N is the 
number of domains in the system being magnetized. The 
same behavior is found in our simulation of a magnetic 
system with nearest neighbor interactions and random- 
ness. A question remains: how large is N for typical 
magnetic tapes? 

For analog storage, we can estimate the number of do- 
mains per cycles. The typical ferromagnetic grain sizes 
found in particulate media magnetic tapes are 0.5 fim in 
length and 0.1 [im in diameter [7| p^p^ |, and the area cov- 
ered by one grain is about 5 x 10^^ cm 2 . The grains used 
in magnetic recording are usually too small to contain a 
domain wall and can therefore be considered as single 
domain particles M. The packing on the tape is usually 
less than 40%. If we assume the percentage to be 35%, 
then the surface grain (domain) density is 7 x 10 s grains 
per cm 2 . The magnetic tapes are typically 1.28 cm wide 
(1/2"), and therefore in 1cm (in length) of tape, there 
are about 9 x 10 s grains. For typical consumer tapes, 
the speed at which the tape is moved is close to hcm/s 
JTj,p3[| . Therefore, the number of grains (domains) per cy- 
cle is 7.5 x 10 7 for a 60 Hz signal, and 2 x 10 5 for a 20 kHz 
signal. Professional tapes have speeds of up to 76 cm/ s 
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and the number of grains N per cycle is 1 x 10 9 and 3 x 10 6 
respectively, for the two frequencies. Therefore, since N 
is large, a V^/V drop in the signal fluctuation is quite sig- 
nificant. (For digital storage, the recording densities are 
as large as 20 million bits per cm 2 |25|,|26|], which for a 
polycrystalline thin film medium gives ~ 5000 grains per 
bit of information, for a 10~ 3 \irri 1 grain size.) As for 
the erasure of the stored information as it is retrieved, 
the added fidelity and linearity should compensate for 
having to rewrite the tape after it is read. 
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